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INTRODUCTION 

This memorandum discusses at some length the five Fortran 
subroutine types. Included are descriptions of the essential 
differences between them, their purpose, definition, and 
proper usage. Special problems associated with the use of 
subprograms are discussed in detail. 

This is essentially a coordination and expansion of 
material gathered from the Fortran manuals published by IBM. 

GENERAL DESCRIPTION 

A subroutine may be considered any pre-defined sequence 
of instructions which performs a specific operation and which 
may be executed repeatedly, at various points within a 
program. 

Fortran provides two kinds of subroutines. One kind is 
the set of functions, of which there are four types. The 
other kind are known as Subroutine subprograms. The functions 
bear a close resemblance in many ways to "functions" in the 
mathematical sense; Subroutine subprograms may also represent 
mathematical functions, but they lack certain resemblances 
and play in fact a much more general role. 

In mathematics, a function is a quantity whose value 
depends upon the values of variables (arguments of the 
function) in a fixed relationship. The function 

f (x) = Ax^ + Bx (1) 

implies a series of operations to be carried out on x, 
regardless of the value of x; however, f (x) itself is. a 
quantity, having a different value for each value of x. Thus, 
one may write 



g(x,y) = ^<^> " ^<y> 
X - y 



(2) 
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illustrating the status of f (x) as quantity , as well as the 
interchangeability of arguments. Lastly, one might define 
a variable z, such that 

z = g(x,y) f(x). (3) 

In Fortran, the various functions have much the same 
status. Equation (1) is in effect the definition of the 
function: this could be translated in Fortran (as one type of 
function) as follows : 

FIRSTF(X) = A*X**2 + B*X. 

Equation (2) illustrates the use of the function, or its 
appearance as a variable, in an arithmetic statement, (which 
is incidentally another function definition). In Fortran, this 
might be written 

SEC0NDF(X,Y) = (FIRSTF(X) - FIRSTF(Y))/ (X-Y) 

while equation (3) would become 

Z = SEC0NDF(X,Y)* FIRSTF(X). 

It becomes clear then that a function in Fortran has 
the following properties : 

1) Its definition is distinct from its use, or call. It 
is called within an arithmetic statement where the function 
value is an operand in the statement. 

2) Hence, the name of a function is the name of its 
value. During execution the instruction sequence that comprises 
the subroutine is carried out, and the result replaces the 
function name in the arithmetic statement. 

3) The arguments given in the definition of a function 
are replaced, when the function is referenced in an arithmetic 
statement, by the arguments required at that time. Hence, 

the arguments in the definition are in effect dummy 
arguments . 

The only major restriction on Fortran functions which 
is not generally applied to mathematical functions is that 
Fortran functions are always single-valued. 

In Fortran, functions, then, are subroutines which 
represent single-valued functions of single or multiple 
arguments, called within an arithmetic statement, where the 
name of the subroutine is the name of the function value. 

A Subroutine subprogram is by nature a more general 
type of subroutine. It may, in fact, be just like a function 
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in its definition, in the status of its arguments^ and in its 
mathematical significance. It may be quite different, however, 
in all these respects, and is always different in the following: 

1) A Subroutine subprogram is called by a separate special 
statement which initiates the execution of the instruction 
sequence which it comprises . 

2) Hence, the name of a Subroutine subprogram corresponds 
not to a value but to the instruction sequence itself. 

In addition, a Subroutine subprogram may produce multiple 
outputs. Its "arguments" may be not only arguments in the 
strict mathematical sense but may be input of any sort from the 
calling program to the Subroutine subprogram; in particular, 
the argument list may include names to be assigned values 
computed by the subroutine. 

FUNCTIONS - DEFINITION 

There are four types of functions in Fortran. The first 
two are used, but not in general defined, by the Fortran 
programmer. The second two are used and usually defined by the 
programmer . 

1. Built-in Functions . These exist in the Fortran 
processor. They differ from all the other subprogram types in 
being open; that is, the instruction sequence is generated by 
Fortran and incorporated into the program, during compilation, 
at each occurrence in the source program of the function name. 

2. Library Functions . These are independent subprograms 
which exist in a program library, either as a card deck or on a 
library tape. They are originally coded in machine language 
according to definite conventions . 

3. Arithmetic Statement Functions . These are defined by 
an arithmetic statement in the calling program. Following their 
definition they are used, or called, by other arithmetic state- 
ments in the same program. 

^' Function Subprograms (also called Fortran Functions), 
These are independent subprograms defined by separate compila- 
tion. The first statement in the definition program has the 
form 

FUNCTI0N name(arg^, arg2, ... arg^) 

which identifies the subprogram as being of function form and 
implies the identification of the value of the function with 
the specified name. 
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SUBROUTINE SUBPROGRAMS - DEFINITION 

These are all of one type. They are independent subprograms 
defined by separate compilation. The first statement of the 
definition program has the form 

SUBROUTINE nameCargj^, arg2, ... arg^^) 

which identifies the subprogram as being of Subroutine sub- 
program form and may imply the identification of one or more of 
the argument names with the value (s) computed by the subroutine. 
The execution of the subroutine is effected by a statement in 
the calling program of the form 

CALL name(arg-j^, arg^^ ... arg^^) 

These five types of subroutines may be differentiated on 
the basis of their form (restrictive or not) and also on the 
basis of the generality of their application (i.e. distinguishing 
functions frequently used by most programs from special- purpose 
routines). One may expect that those most commonly used should 
be most readily available, or simplest to use. 

1. The built-in functions are a special case, being open 
routines. In fact, they look like subroutines from the vantage 
point of the source program but not from that of the object 
program. Since the generated coding is reproduced in the 
object program each time the function is referenced, such .. . ; ' 
routines must be limited to functions which are simple in terms 
of machine- language coding. There is a decrease in execution time 
in the object program, though at the expense of memory space. 
Since these routines are implicit in the Fortran Processor 
(occupying valuable space and time in it) they must be restricted 
to functions which are of the most general utility and which 
become an extension of the Fortran language itself, rather than 
being germane to its applications. For such basic functions, the 
open type combines the simplicity of subroutine usage in the 
source program and the efficiency of machine- language straight 
line coding in the object program. 

2. Library functions, as the- name implies, are of suf- 
ficiently general utility to be permanently maintained by the 
computer installation. The most commonly used of these are 
maintained on a special tape which is automatically available 
to every program running in the system. Others are maintained 
as binary decks also readily available to the user. Functions 
made available in this form are likely to include those most 
often needed by the user of a particular installation -- their 
generality, in other words, is peculiar to the principal appli- 
cations of one computer. 

Library functions, since they are pre-written and used 
frequently, are coded in machine language mainly in the interests 
of economy. They are subprograms with a unique convention for 
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linkage which demands machine- language coding, 

3. Arithmetic statement functions are the first and 
simplest of the subroutines written by the user. They are 
restricted to those functions which may be expressed in a single 
Fortran arithmetic statement. They satisfy, however, the basic 
requirements of a subroutine — i.e., repeated availability 
without redefinition, and interchangeability of arguments. 
The limitations on their definition are counterbalanced by a 
distinct advantage over subprograms: since they are internal 
to the calling program, all the variables of the calling 
program are available to the subroutine as arguments and 
variables, without the problems of transmission introduced by 
the use of other subprograms. 

^« Functio n subpros;rams are used for the computation of 
functions which may not be defined by a single arithmetic state- 
ment. They may, therefore, be mathematically as complicated as 
necessary. (We may note here that at one end of the scale are 
subroutine types which are highly restrictive but which compute 
functions of the most common type; on the other end are sub- 
routine types of the most generalized form, used for more special 
applications . ) 

Function subprograms are principally used for mathematical 
operations. This is implicit in the major restriction placed on 
them, that they compute a single quantity which takes its place 
in the arithmetic statement which referenced the subroutine. It 
is possible to extend the application of these programs so that, 
for example, they may compute many other quantities as well, or 
may perform additional non-arithmetic operations. Their most 
usual application, however, and the one for which they are. 
intended, is the computation of a single-valued arithmetic 
f unc tion . 

5' Subrouti ne subprograms . finally, are the most general 
form of subroutines available. They can be used to do anything 
a Function subprogram does, but are even less restrictive. 
Probably their most fitting applications are for the computation 
of multi-valued functions (although this can also be done in- 
directly by means of a Function subprogram) and especially for 
the execution of non-arithmetic operations (e.g. input and 
output). 

FUNCTIONS - NAMING AND CALLING 

1. Built-in, Library and Arithmetic Statement functions are 
all named and called in the same manner. 

The name consists of 4 to 7* alphanumeric characters with 



Double precision and complex Library and Built-in functions: 
4 to o characters . " 
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these restrictions: the last character must be F; the first 
character must be alphabetic^ and is X if and only if the value 
of the function is of fixed point mode. 

N.B. The mode of a function and of its arguments is determined 
by the function definition. Hence, a change in the mode configu- 
ration would require a new definition and a new name.** 

Examples ; 

Built-in 

I = (Mi-N/2)*XM0DF(M,N) 

X = DELTAX + MAX1F(YA,YB,YC,10.0) 

Library 

B = A + C0SF(A) 

M = XDETSMF(NN,N,AMATR,D) 

Arithmetic Statement 

(Definition) TESTF(ARG1,ARG2) = (A/ARG2)*(ARG1-ARG2) 



(Use) A = 3.0*DELTA 

GAMMA = ALFA*TESTF (ALFA, BETA) 



2. Function Subprograms . The name consists of 1 to 6 
alphanumeric characters with these restrictions: the last 
character must not be F unless the entire name is less than 4 
characters; the first character must be alphabetic- and .must 
govern the. mode of the function according tc the rules for « 
naming variables. This name is used: to call the .function in the 
same ..manner as the ■ other f unc tions . No E . is attached to the 
function name when it is used in a statement. 

Examples ; 

(Definition) FUNCTI0N FUNK(A,B,C) 



FUNK = P+Q*B*C 
RETURN 



''* There are no fixed-point double precision or complex function 
names. The names of double- precision and complex Lxbrary and 



Built-in functions are listed with a prefix D or V. This prefix 
must noiE' be used m function references. 
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(Use) . 



END 

A = FUNK(DELTA,X,25.) + ALPHA 



SUBROUTINE SUBPROGRAMS - NAMING AND CALLING 

The name consists of 1 to 6 alphanumeric characters with 
these restrictions: the first character is alphabetic; the last 
character must not be F unless the entire name is less than 4 
characters. Since the name of the subroutine is not the name of 
a value, the name does not affect the mode of any of its values. 

The names of values computed by the subroutine may be 
designated in the argument list. These govern the mode of the 
values they represent. There must be correspondence in the mode 
of each argument in the subroutine definition with the respec- 
tive arguments in every call for the subroutine. 

EXAMPLES : 

(Definition) SUBROUTINE L0VE (ABLE, BAKER, NAN) . 



(Use) 



NAN = ... 

RETURN 

END 

X = X + DELX 

CALL L0VE(X,Y,M) 

Y = Y**M 



DOUBLE-PRECISION AND COMPLEX FUNCTIONS AND SUBROUTINES 

An arithmetic statement referring to a double -precis ion or 
compleK function is prefixed with a D or I in column 1. AD 
or I prefixing a CALL statement signifies that all floating- 
point arguments in the statement are double-precision or complex, 
respectively. 

A double- precis ion or complex Arithmetic Statement function 
definition must of course have D or I in column 1. 



EXAMPLE : 



D AREAF(R) = SQRTF (PI*R**2) 



USE OF FUNCTION AND SUBROUTINE NAMES AS ARGUMENTS 



The argument list of a Function subprogram or a Subroutine 
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subprogram may include the names of Library functions, other 
Function, subprograms, or other Subroutine subprograms (but not 
Built-in or Arithmetic Statement functions). As with any other 
argument, a dummy function or subroutine name appears in the 
argument list of the subroutine definition. The following rules 
apply: 

1. The name of a Function subprogram. Subroutine sub- 
program, or Library function used in an argument list must 
appear also in a pseudo- statement, in the calling program, of 
the form 

F NAMEl, NAME2, .., NAMEn (F in column 1). 

2. In the case of a Library function the terminal F is 
dropped in the argument list and also in the F pseudo-statement. 

3. In the case of a double-precision or complex Library 
function, the name appearing on the F card is prefixed with a D 
or I, 

EXAMPLES ; 

SUBR0UTINE PRIME (DUMMY, X, A) 



A = DUMMyF(X) 

• « • 

This would allow the selection of any library function when 
calling PRIME, as 

CALL PRIME (SIN, ALPHA, Q) 

• « • 

CALL PRIME (C0S, ALPHA, P) 

• a • 

F SIN, C0S 
Again, 

FUNCTI0N PETER(DUMMY,X,A) 

• • • 

■ • • • 

PETER = DUMMY (Y) 



allows the calling program to specify various Function sub- 
programs when calling PETER, as 
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X = Y + PETER(FIRST,ARG,Ai) 
A = PETER(SEC0ND,ARG,A2) 

* • « 

F FIRST, SEC0ND 
Finally, 

SUBR0UTINE PAUL ( DUMMY, X, A). 

• • * 

CALL DUMMY (X) 



allows the calling program to specify subroutine subprograms ■. :. 
when calling PAUL, as 

CALL PAUL(SUBR,EX,EA), 

• • « 

F SUBR 

Note that the definitions (though not the argument list) of 
Arithmetic Statement functions, as well as of Function sub- 
programs and Subroutine subprograms, may involve other functions. 
An Arithmetic Statement function definition may involve other 
Arithmetic Statement functions which have been previously : 
defined in the same program. 

RULES FOR DIMENSIONING SUBROUTINE ARGUMENTS 

The arguments in a SUBR0UTINE or FUKCTI0N statement (the 
du3imiy arguments of the subprogram definition) may not be 
subscripted. No subscripted variables may appear in the state- 
ment defining an arithmetic statement function. An array may be 
an argument; the array name (non-subscripted) appearing in a 
CALL statement or function subprogram reference must of course 
be dimensioned in the calling subprogram; the corresponding 
dummy argument of the subprogram definition must be dimensioned 
in the subprogram and the dimensions must correspond in both 
calling and called programs. 

In fact, the dimensions of the dimimy array in a subprogram 
need not be exactly the same as those of the corresponding 
array in the calling program. Variation is permitted as follows : 

The dimension of a one -dimensional array may vary; 

The second dimension of a two-dimensional array may 
vary ; 

The third dimension of a three-dimensional array may 
vary; 
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in other words, for an n-dimensional array, the first (n-1) 
dimensions must correspond exactly between subprograms. 

A single element of an array may be used in the argument 
list of a calling statement. It is, of course, subscripted, and 
corresponds to a single non-subscripted dummy variable, as : 



Calling : 
CALL SUB(A,B(I,J)) 
On the other hand. 

Calling : 
DIMENSION B(15, 10) 



Called: 
SUBROUTINE SUB(X,Y) 

Called: 
SUBR0UTINE SUB(X,Y) 



CALL SUB(A,B) 
transmits the whole array B, 



DIMENSI0N Y(15, 10) 



An exception to the above occurs in the case of double- 
precision or complex variables. Given a double-precision or 
complex array (say ALPHA, dimensioned (n,n)) where it is 
desired to use a single element ALPHA(I,J) as an argument to a 
subprogram, then the corresponding dummy argument in the sub- 
program must be given like dimensions (n,n) in a double- .. 
precision or complex DIMENSI0N statement. Otherwise ALPHA (I, J) 
must be set equal to a single variable for use in the argument 
list. For example, the following illustrates- both cases: 



Calling : 

D DIMENSI0N A(10),B(5,5) 

D C = A(5) 

D CALL NAME (C,B(1,3)) 



Called : 

SUBROUTINE NAME(X,Y) 
D DIMENSI0N Y(5,5) 
D Y = X + (2.,1.) 

RETURN 



Dimensioning Y in the subroutine enables the FORTRAN processor 
properly to compute the address of the imaginary or least signi- 
ficant element of the pair represented by B(l, 3). Reference is 
then made in the subprogram to the variable Y, which is in fact 
the first element, i.e. Y(l,l), of an implicit array. Use of 
subscripts (other than (1, 1)) with Y should be employed only 
with great care and a thorough understanding of the ordering of 
arrays in storage. 



USE OF C0MM0N FOR SUBPROGRAM LINKAGE 

Variables may be transmitted to a function or subroutine 
subprogram by the use of C0MM0N storage. Consider the arithmetic 
statement function 
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FIRSTF (X,Y) « (A*X+B)/Y . 

The arguments X and Y are dummy argximents and will be replaced 
when the function is called; A and B, however, are variables 
x^hich always remain the same (though not necessarily of the same 
value). This sort of arrangement is possible because the 
arithmetic statement function is, though closed, contained with- 
in the calling program. 

In the case of a function or subroutine subprogram, however, 
where the subroutine is an independent program and all the 
programs are relocatable, variables like A, B which we may call 

function constants", must be transmitted from the calling to the 
called subprogram. One could write 

CALL SUBR (X,Y,A,B) 
CALL SUBR (P,Q,A,B) 
etc., 

but it is also possible to use C0MM0N storage for the purpose. 

The statement 

C0MM0N A, B, . . . . 

assigns the given variables to specific non-relocatable loca- 
tions in upper core. The first variable is assigned to 77461g, 
the second to 77460g, and the rest in order proceeding down- 
ward. Successive C^MM0N statements x^ithin a program will 
continue to store variables below those already encountered. 
Each subprogram will again assign the variables listed in C0MM0N 
statements to upper core starting at 77461g. Thus one may write: 

Calling: Called: 

DIMENSI0N T(5,5) SUBR0UTINE SUBR<D) 
C0MM0N A, B, T, C DIMENSI0N DUMMY (5, 5) 

C(&MM0N P, Q, DUMMY, R 
CALL SUBR(X) 

Because of the correspondence in the order of variables in the 
two C0MM0N statements, the pairs A and P, B and Q, and T and 
DUMMY will be identical. 

For most purposes, the simplest and safest way to guarantee 
consistency among subprograms is to make C0MM0N statements, and 
DIMENSION statements for those variables which appear in C0MM0N 
statements, identical in all subprograms and the main program. 
This may easily be done and there is, of course, no need, 
generally, to change variable names between subprograms . (The 
use in those examples of different names for arguments in sub- 
program definitions is primarily to stress their status as 
dtmmy variables.) As will be seen from what follows, departures 
from this procedure should be made with caution. 
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In the example above, should it be desired to call another 
subprogram which involves the variable C but not the others, it 
would be necessary to include in a C0!yiM0N statement in the sub- 
program, dummy variables of proper dimension to assign C0MM0N 
storage properly, thus: 



SUBROUTINE NAME(ARG) 
DIMENSI0N DUMMY(5,5) 
C0MM0N P, Q, DUMMY, R 



(1) 



This is necessary to force correspondence between all CifiMM0N 
assignments. 

Additional care must be taken with double-precision or 
complex variables assigned to C0MM0N storage, where several sub- 
programs refer to C0MM0N storage. Consider the example of a 
main program where A and B are complex variables and T is a 
complex array: 

I DIMENSION T(2,4) 

C0MM0N A,B,T 
I ... 

I X = A**2+B**2 
I CALL ABLE(X) 



Since complex variables consist of two parts, they are stored 
like arrays. In this example, C0MM0N storage would look like 
this : 



77461 
60 
57 
56 
55 
54 
53 



A (real part) 

A (imaginary part) 

B (real part) 

B (imaginary part) 

T(l, 1) (real part) 

T(2,l) (real part) 

T(l, 2) 



77445 : T(l, 1) (imaginary part) 
etc. 

The appearance of A and B in complex arithmetic statements, and 
of T in a complex dimension statement, guarantees that double 
storage will be assigned. Suppose now that the subprogram 
Able is to operate on the array T, but does not involve A or B. 
If we write 
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SUBR0UTINE ABLE(D) 
I DIMENSION DUMMY (2, 4) 
C0MM0N P, Qj DUMMY 



then P and Q are for the purpose of "spacing down" to the loca- 
tion of the array in upper core. Since we have assumed, however, 
that P and Q never appear again in ABLE, there is within. this 
independent subprogram no indication that P and Q are complex 
variables. Hence C0MM0N storage will look like this: 



77461 
77460 
77457 
77456 



P 

Q 

DUMMY (1,1) (real part) 

DUMMY (2,1) (real part) 



etc, 



and the desired correspondence between the arrays T and DUMMY is 
not accomplished. 

The simplest way to circumvent this problem is to include 
all complex (double-precision) variables assigned to C0MM0N in a 
complex (double-precision) DIMENSI0N statement. Thus, in the 
example above, if we write 

SUBR0UTINE ABLE(D) 
I DIMENSI0K DUMMY(2,4), P(l), Q(l) 
C0MM0N P,Q, DUMMY 

The dimensions of P and Q are not really changed but their status 
as complex variables is established, and each will be assigned 
two locations , 

Clearly, the use of C0MM0N assignments may be applicable 
not only to what we have called function constants, but to 
arguments as well. Arguments which are transmitted to a sub- 
program by being placed in C0MM0N storage are implicit 
arguments, as opposed to those explicitly mentioned in the sub- 
program argument list. Again, it is only necessary to create 
correspondence in upper core between each argument and its dummy 
correspondent in the subprogram. 

It must be noted that a Function subprogram must have at 
least one argument. Thus, implicit arguments may be used for all 
argtraients of a Subroutine subprogram, but at least one explicit 
argument (a dummy, if so desired) must be given for a Function 
subprogram. 



-14- 

RELATIOM BETWEEN CUmM AMD EQUIVALENCE STATEMENTS 

The normal order of Cj5MM0N storage assignment, described 
above, is changed in the case V7here variables appear in both 
C^MM0W and EQUIVALENCE statements. Any CJ&MM0K variables 
appearing in EQUIVALENCE statements will be assigned upper 
core values above those not in EQUIVALENCE statements 
regardless of the order of the C0MM0N statements. For example, 

C0MM0N A, B, C, D 

EQUIVALENCE (C,G), (E,B) 

will cause the following assignments : 

77461g C and G 

60 B and E 

57 A 

56 D 

Note that the variables are assigned in the order in which 
they appear in the EQUIVALENCE, and then in the order of the 
C0iMI'I0N statement, (Note also that a variable made equivalent 
to a C0M0N variable is assigned to C0MM0N storage.) 

Consequently, if equivalences are established involving 
variables in C0MM0N storage, similar equivalences should be 
sat up in all subprograms which refer to C0MM0W in order to 
force proper correspondence between subprograms. 

In the case of double-precision or coraplesc variables, 
a D or I prefixing a C0MM0N or EQUIVALENCE statement is ignored, 
but the double nature of the variables is preserved when they 
are so defined in other statements vjithin the program. Generally, 
single-precision variables should not be made equivalent to 
double-precision or complex variables, nor should they be 
given the same location in C0MM0N storage. 

TRANSMISSION OF VARIAB L E SUBSCRIPTS BETOSEN SUBPROGRAMS 

A variable subscript has a somewhat ambiguous status 
because of the way FORTRAN performs indexing. If a D0 loop 
(or an input list) has an index I which appears only as index 
and subscript, there is no storage location set aside for I: 
the indexing and the subscript reference are taken care of in 
another vjay. If I also appears as a variable, either inside 
or outside of the loop, then there is a storage location 
assigned to I, and this location Is referred to when necessary. 
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As a result, the storage location "I" may sometimes have a 
value different from the indexing value of I. For example: 

I = N/2 
A(I) = X 
D0 1 I = 1,M 
B(I) = X*A(I) 
1 CjONTINUE 

When the loop is completed, i.e., when the index I is equal 
to m-1, the storage location *'l" will still contain the value 
N/2, and any reference to I as a variable or as a subscript 
will use this value until a nev? value is assigned or a new 
loop is initiated. Suppose we continue the above program, 
following statement 1 by: 

NU = 1**2 (equivalent to NU = (N/2)**2) 
D0 2 I = 1,N 

• • • 

2 L =2*1-1 (equivalent to L = 1,3, ... 2*N-1) 

Following statement 2, the storage location of I (i.e., I in 
the role of variable ) will be equal to NH-1, the current value 
of the index I. This is because I is used as a variable within 
the loop. 

Because of the ambiguous nature of subscripts, care must 
be exercised in transmitting them between subprograms. (N.B.: 
A variable subscript x^hich is not being used as an index is 
called, in the FORTRAN manual, a ''relative constant". For 
further details about variable subscripts, see Fortran 
Reference Manual, pp. 82-35.) 

1, In an Argument List 

Transmission of variable subscripts via an argument list 
is straightforward. 

I «= N/2 SUBRJJUTINE SUBK.(A,N,Y) 

CALL SUBR(X,I,B) DB1ENSI0N Q(10) 

Y = Q(N) 

The value of I will be transmitted properly. Other examples 
of correct linkage are: 

D0 1 I « 1,10 
CALL SUBR(X,I,B) 
1 ... 

where although I is not explicitly used as a variable its 
appearance in the argument list allov7S for proper transmission. 
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and, 

CALL SUPER (I,Y) SUBR0UTINE SUPER (K,A) 

K = XFIXF(A**2-X**2) 
RETURN 

where I will be properly transmitted from the calling to the 
called program and back again. Note^ however, that if the 
CALL statement in this example were within the range of a D0 
with index I, the change in I effected by SUPER would not hold, 
and I would be reset, upon return, to the proper indexing value. 

It should be observed that in this last example, the 
argument K is used both as an input and an output parameter. 
This is in general a dangerous procedure; any change in the 
value of input variables may cause difficulties, because it 
requires considerable understanding of FORTRAN operation to be 
sure of exactly how this change will be effected. 

2. Transmission of Variable Subscripts via C(&MM0N Statements 

Consider the following example: 

C0MM0N K (1) SUBR0UTINE ABC(X) (6) 

K = K + 1 (2) 

A(K) = B(K)*^V2 (3) C0MM0N N (7) 

CALL ABC(A(K)) (4) 

TERM = B(K) (5) N = N + 2 (8) 

When ABC is called, K has its proper value in C0MM0N 
storage and thus N has the same value. Statement (8) changes 
the value of K in C0MM0N storage. At execution of statement (5) 
however, the value of K in C0iMM0N is not used to find the 
appropriate term in the array B; because of the way FORTRAN is 
set up the subscript K is determined by a different storage 
method, and on return from a subprogram which recalculates the 
variable K, the subscript K will not have been changed. 

The statement 

K = K (4a), 

inserted between (4) and (5) in the example, will cause the 
subscript K to be reset by the value of the variable in C0MM0N. 

Now consider this example. 
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C0MM0N K (1) SUBR0UTINE ABC(X) (4) 

C0MM0N N (5) 

D0 1 K = 1,J (2) 
1 CALL ABC(A(K)) (3) 

Throughout the loop on K, the subscript K will be 
advanced in accordance with the D0 specification; the value of 
K in C0MM0N storage, however, as transmitted to the subprogram 
ABC, will remain at whatever value it had before entering the 
loop, (If statement (2) is the first appearance of K following 
statement (1), moreover, K in C0MM0N storage will have had no 
value assigned to it at all.) A statement 

K = K (2a) 

inserted within the range of the D1&, or any appropriate state- 
ment with K on the left of an = sign, will reset K in C0MM0N 
storage (i.e., K in its role as variable ) in accordance with 
the values of the index. 



